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The production of 77')-mesons during the decay of a CP-odd phase is studied within an 
evolution operator approach. We derive a quantum kinetic equation starting from the Witten- 
DiVecchia-Veneziano Lagrangian for pseudoscalar mesons containing a Ua{1) symmetry breaking 
term. The non-linear vacuum mean field for the flavour singlet pseudoscalar meson is treated as 
a classical, self-interacting background fleld with fluctuations assumed to be small. The numerical 
solution provides the time evolution of momentum distribution function of produced rj'- mesons after 
a quench at the deconflnement phase transition. We show that the time evolution of the momentum 
distribution of the produced mesons depend strongly on the shape of the effective potential at the 
end of the quench, exhibiting either parametric or tachyonic resonances. Quantum statistical effects 
are essential and lead to a pronounced Bose enhancement of the low momentum states. 
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I. INTRODUCTION 

Construction of the Relativistic Heavy Ion Collider 
(RHIC) at the Brookhaven National Laboratory is com- 
pleted and it is designed to initiate energy densities suffi- 
cient to produce a quark gluon plasma (QGP) Such a 
strongly correlated state of matter has a life time smaller 
than 1/m/c and due to rapid collisions the plasma ther- 
malizes, and at critical values of temperature and den- 
sity the quarks and gluons form hadronic bound states: 
a process driven by confinement and chiral symmetry 
breaking. Many aspects of the plasma's production and 
evolution are characterised by non- linear dynamics. The 
hadronisation process itself as well as critical phenom- 
ena in the vicinity of the phase boundary require a study 
with non-equilibrium techniques. 

One challenging example of a far from equilibrium pro- 
cess is spontaneous particle creation in a strong back- 
ground electric field, i.e. the Schwinger mechanism [Q. 
However, pair creation in QED has never been observed 
directly although planned new facilities such as a X-ray 
free electron laser (XFEL) [^jjij will allow to reach the 
region of required critical field strengths. Therefore this 
non-perturbative effect was mainly studied for applica- 
tions providing strong enough fields ranging from black 
hole quantum evaporation |^ to particle production in 
the early universe |^ and in ultra-relativistic heavy-ion 
collisions 

An unsolved problem of conceptual and practical in- 
terest is the precise connection between field theory and 



kinetic theory. Recently a link between the mean field 
approach of vacuum pair creation in a spatially homo- 
geneous Abelian background field and a kinetic formula- 
tion was established in Q . The resulting source term for 
spontaneous pair creation is non-Markovian and retains 
quantum statistical effects [^,0. In many approaches 
the background field is treated as a time dependent clas- 
sical field with feedback incorporated via Maxwell's equa- 
tion, e.g. 1 11 15). In these approaches the production of 
fermion/gluon pairs was employed to describe the for- 
mation of a quark-gluon plasma. Herein we focus on the 
production of bosonic particles in hot hadronic matter in 
QCD. 

Lattice calculations, e.g. [pj||, as well as QCD Green 
function approaches, e.g. |15| , |17| ], indicate that the de- 
confinement and chiral phase transitions are coincident 
]l5| , p^ [2C|] . At present it is an open question whether 
the restoration of the Ua(1) symmetry which is broken 
in the QCD vacuum sets in already at the deconflne- 
ment transition temperature or above. In addition par- 
ity maybe spontaneously broken which is connected with 
a non- vanishing QCD angle [|l). The CP odd phase 
is of particular interest since it may have experimental 
signatures such as an enhanced production of 77 and 77' 
mesons which can contribute via their decays to 

the low mass dilepton enhancement. They can decay via 
CP violating processes such as 77 — > vr^vr". 

Herein we study the production of rj'- particles dur- 
ing the decay of the CP-odd phase. Complementary to 
p3 where the decay rate of metastable states was esti- 
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mated, we study the full time evolution of the momentum 
distribution function using a kinetic description based 
on an effective Lagrangian. We start from the Witten- 
DiVecchia-Veneziano model |25|] , however, different ap- 
proaches can be applied, e.g. ]2(|. 

In this article, the external background field concept 
is replaced by a potential yielding self-interaction and 
non-linearity. This potential dominates the solution of 
the quantum kinetic equation which is derived using an 
evolution operator approach. The introduced technique 
to link an effective Lagrangian and kinetic theory is not 
restricted to the discussed model calculation of 77' produc- 
tion. It's application is general in quantum field theory. 

The article is organised as follows. In Section II we 
introduce the model Lagrangian and identify the self- 
interaction parts. In Section III we perform the quan- 
tization of the evolution operator used in Section IV to 
derive a quantum kinetic equation. In Sections V and VI 
we discuss the decay of the CP odd phase in view of our 
numerical results. 



/i and a are functions of T and hence the potential cor- 
responding to (H) has modified properties close to the 
deconfinement phase transition. 

From we obtain the following Klein-Gordon type 
equation of motion for the field rj{x^t): 



where TOq 



(□-|-m^)ry = J,, 
^ The nonlinear current 



(3) 



(4) 



contains orders rf and higher and is related to the self- 
interaction of the field 77. Note that the linear term of 
the total current J — — /i^?/ + Js is contained in the mass 
squared term of the lelf-hand side of Eq. 

The total Hamiltonian density, H. = Tio + 'Hs , is given 

by 
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II. THE EFFECTIVE LAGRANGIAN 
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We start from the effective Lagrangian of the Witten- 
DiVecchia-Veneziano model |^5| 

-Ce// = ^ {tr(^^JJ^^JJ^) + tr{MU + Af (7+) 



[6*- -ir(lnC/-lnC/+)l2 



(1) 



where TYq involves only the free field part with the mass 
Too; Tis includes self interaction starting at orders 77^ and 
TT is the momentum canonically conjugate to 77: 



7r(f, t) = 77(x, i), 



(6) 



where the overdot denotes the derivative with respect to 
time. 



which describes the low-energy dynamics of the nonet 
of the pseudoscalar mesons |^ in the large A*"c-limit of 
QCD. The meson fields are described by the N f x 7V/- 
matrix \J in Eq. (Q). Explicit chiral symmetry breaking 
is realized by the current quark mass matrix M with 
the diagonal elements related to tt and K meson masses. 
With the parametrization U = eKp(i(f>/ f^^) , the matrix 
(j) representing the singlet and the octet meson fields 
yields the pseudoscalar nonet. The last term in the ef- 
fective Lagrangian is related to the ?7A(l)-anomaly: the 
singlet is massive also in chiral limit. The parameter 
a = 2NfXYM/f^ contains the topological susceptibil- 
ity, Xym- Herein we focus on the singlet state which 
is the main component for 77' and obtain the following 
Lagrangian: 
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(2) 



III. EVOLUTION OPERATOR APPROACH 

We introduce the in-field, 77^,1 (af, as a solution of 
Eq. in absence of sources and quantise it according 
to the standard canonical procedure (see Appendix A). 
The original self-interacting field is connected with the 
in-field by the unitary transformation: 



where 



U{t) = Texp{^i f dt'Hi"{t')} 

J —00 



(7) 



(8) 



is the time evolution operator with the self-interaction 
Hamiltonian written in terms of the in-held operators 



In Eq. 



/ = i/I/tt, where f^^ — 92 MeV is the semi- 



leptonic pion decay constant, = ^{m'l + 2mj^) is a pa- 
rameter depending on tt- and if-meson masses. For zero 
temperature T = 0, a to^ + to^, - 2to|. ~ 0.726 GeV^ 

and /x^ ~ 0.171 GeV . In response to non-zero tempera- 
ture and density mesons have an effective mass, e.g. [p8[: 



^The model is defined in a finite volume: V — L^, — L/2 < 
Xi < L/2, i = 1,2,3. The continuum limit is jr'^j: => 

J (2^- 
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i7r = J d^xnsiv^V^n■,7^^^r,„). (9) 

In the limit t — > —oo we have U{t) I, so that 

hm ri{x,t) ^ ri,n{x,t). (10) 

t — ^ — oo 

The exact meaning of (|l^) depends on details of the 
current which in our model is determined by self- 
interaction taking place at all times. Hence Eq. ( |l0| ) 
is a priori difficult to justify. We assume an adiabatic 
vanishing of the interaction for i — > — oo. 

The field ri{x,t) is given by the space- homogeneous 
mean value 4'{t) — {ri{x,t)) and fluctuations x 



T]{x,t) ^ (f>{t) + X{x,t) 



(11) 



with {x{x,t)) — 0. Assuming that X^ f ■• quantum fluc- 
tuations can be treated perturbatively. Herein we restrict 
ourselves to zeroth '^"^ and flrst order. Substituting 
Eq. (|ll|) into Eq. (|) yields 



{U+ml)x + (t> + ml4>^ 



S 1 



where 



1 — cos 



(I) 



The zeroth order of the current is given by 



(12) 



(13) 



(14) 



Taking the mean value 
mean field equation 



of Eq. (f2|) , yields the vacuum 



ffi sin 



= 0. 



(15) 



Eq. (jlj) in concert with Eq. (|l2|) provides the equation 
of motion for the quantum fluctuations 



(□+mg)x-M'(l-cos(|) 



(16) 



The right-hand side of this equation vanishes in the in- 
limit. Rewriting ( p^ ) for the Fourier components xik,t), 
we obtain a Mathieu type equation |29,p(l| 



where 



Xik,t)+iulit)x{k,t)^0, 



.^(t)^(c.°)^-^2(l-cos(|)) 



(17) 



(18) 



and iOkit) is the time-dependent frequency of the fluctu- 
ations with limt^_oo 0Jk{t) = uj^ = y^W + m^. 

For a > fj.^ , the frequency squared is positive for all 
momentum modes and at all times. However, if a < /i^, 



tj^(t) can be negative for modes below a critical momen- 
tum kc indicating a tachyonic regime. 

It is important to observe that Eqs. ( [l5| ) and (|6|) are 
coupled. Although the fluctuations do not react on the 
vacuum mean field, the latter modifies the equation for 
fiuctuations via a time dependent frequency. 

The self-interaction Hamiltonian density correspond- 
ing to the equations (|l^) and ( |l6|) is quadratic in x, 



+ V 

1 2 



»(f)-(|) 

COS (^) - 



(19) 



and also vanishes when t — > — oo. Hence in the approxi- 
mation of preserving quantum fluctuations in the vacuum 
mean field but neglecting the feedback, the adiabatic hy- 
pothesis of vanishing interactions for t — > — oo discussed 
with Eq. ( p^ is justified. 

For the Fourier components of the fluctuations, we 
write the ansatz analogous to (|A4|), 



X{k, t) = VAt)a{k, t) + n(i)at (_fc, t), 



where 



r£(^) 



exp{-iefc(wfc,i)}. 



(20) 



(21) 



and Ofc(wfc,i) is a phase which in the in-limit takes 
the form Lo^t. In the same limit, ^^{t) -'^^(^)' 
while the time-dependent operators a{k,t), a^{k,t) with 
limt_,_oo a{k, t) = ainik) and limt^_oo a'l'(fc, t) = aj„(fc). 

In the case when the fiuctuations and the frequency ojk 
vary adiabatically slowly in time, the dynamical phase 0^ 
can be chosen as 



ujk{t')dt'. 



(22) 



The relations between the Fourier components r]{k, t) and 
x{k, t) and the corresponding conjugate momenta are 
given by 

il{k,t)=x{k,t) + 5j:^^Vvm. (23) 
TT{k,t)=^^{k,t) + 5^^^^^{t). (24) 

The Fourier components of the operator tt^ are 
n^{k, t) = -iuok{t) [rg(t)a(-fc, t) - rl{t)a\k, t)] (25) 



and in the limit t ^ —oo this ansatz reduces to (A5). 

Using Eqs. ( |20| ) and (^), we obtain the following re- 
lations between a{k,t), a^ik^t) and the in-opcrators: 
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a{k, t) 



a\k, t) 



——{U'\t)[rj,n{k,t) + — ^,„(-fc,i) U{t)) 



\ UJQ J 



(26) 



— i--{C/-l(t) Ln{-k, t) - —n,n{k, t)] U{t) 



(27) 



where ojq = LUk=a- It is easy to verify that the operators 
a{k,t), a^(fc,f) fulfiU the sam e commutation relations as 
the in-operators, Eqs. (|A6|). Hence the transformation 
defined by the evolution operator is canonical. The (f>- 
dependent terms in Eqs. (p^)- ( p7| ) act like counter terms 
which cancel the vacuum mean field contribution of the 
previous terms, so that ( p6| ) and ( |27| ) do not depend on 
(j) explicitly. 



not change in time. One considers excitations with re- 
spect to this vacuum and interprets an increase in the 
occupation number density as particle production. The 
vacuum state itself is "empty", i.e. without particles. 

In our model, the background field (/)(<:) is periodic 
in time, i.e. in addition to the quantum fluctuations 
around (p we have oscillations of (p itself. Therefore the 
vacuum restructures itself at each moment in time and 
consequently the occupation number density has to be 
redefined as well since it is assumed to be zero only for 
the vacuum state. In the time evolution of the densities 
Af±{k, t), it is therefore necessary to separate the contri- 
bution of the real particle production from the one related 
to the vacuum state redefinition. This is achieved by us- 
ing the expansion ( |l^ in the evolution operator U{t). 

We consider first the time evolution of A/'_(fc,t). Tak- 
ing the time derivative of J\f^{k,t) and taking into ac- 
count the relation ill — Hl^U we find: 



IV. KINETIC EQUATION 

The number of particles of a given state characterized 
by the momentum k at time t is given by 



7V(/c',t) = {Q\a\k,t)a{k,t)\Q). 



(28) 



In the limit t —oo, Af^k, t) tends of course towards the 
occupation number density of the in-field: 



U{k,t) N{k) = (0|aL(fc)a»n(fc)|0) 



(29) 



Substituting (^6|) and ( |27D into ( |28| ) and introducing the 
instantaneous states U\0) = \U), {0\U~^ = {U\, the par- 
ticle number can be written as 



AA(fc, t) 



1 



{U\riUk,t)mn{k,t) + ^„,(/c,i)7rJ„(fc,i)|;7) 



-{U\r,l{k,t)^Uk,t) ~ nUk,t)mn{k,t)\U) 



(30) 



The number of particles of momentum k is not equal to 
that of momentum (— fc) for all times t. Therefore it is 
convenient to introduce 

M±{k,t) = ]^(U{k,t)±U{~k,t)), (31) 

where M-\-{k, t) is the particle number averaged over the 
directions k and {—k), while J\f^{k,t) measures the de- 
gree of asymmetry. At fixed volume, the occupation num- 
ber densities can change in time for two reasons, either 
with the change of the number of particles or with the 
change of the vacuum state. The presence of the back- 
ground field leads to a restructuring of the vacuum state. 
Note that in the case when the background is a con- 
stant classical field, the definition of the vacuum does 



J<f-{k,t) = 



(32) 



Using the expansion (19), the commutator in (^2|) is read- 
ily calculated: 

U.{k,t) - -^Im J d'x{0\e''^^x{k,t)Ji'H^,m- (33) 



The time evolution of the density A/"- (fc, t) is determined 
by the self-interaction of the field r]{x, t). To get an exact 
formula valid in all orders of perturbations in (x//) it is 
sufficient to replace ji^^ by Jg in Eq. (|33|). 

In the chosen approximation of small quantum fiuctu- 
ations, the current Js is considered in first order in x- 
Using Eq. (^3|), the integral in Eq. ( |3^ ) turns out to be 
real and one obtains Af-{k, t) = 0. In first order in x the 
number density Af-{k,t) is therefore conserved. 

Taking the time derivative oi Af+{k,t) we obtain the 
evolution equation 



J^+{k,t) = — Re[c(fc,i)e-2'e* 
Wfc L 

1 



i 

2^1 



{Lot~{JlY)lm C{k,t)e 



0^0 



HI 



.(fc,t)4„(fc,i)l \U), 



(34) 



where we have defined the time-dependent pair correla- 
tion function C{k,t) = {0\a{~k,t)a(k,t)\0) . Calculating 
the commutator in the expression for A/+(fc, t), we find 



2ujk 
1 1 
V wfc 



HI 



,(fc,i)7rl„(fc,<) \U) = 



Re / d^x{0\e'''^T:Uk,t)4^^x,t)\0) 



ujo 



(35) 
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With Eq. (|l3|), this last expression becomes 



— {U\[Hl\7r.4it)nl{k,t)\_\U) 



2a;, 
1 



Im 



C(k,t)e 



-2i0. 



(36) 



so that in the smaU quantum fluctuations approximation 



JV+{k,t) = —Re 



C(fc,i)e 



(37) 



In the same approximation, the pair correlation function 
C{k,t) obeys the equation 



c{k,t)^2i{ek--L0k)cik,t) 

|^(l + 2AA+(fc,i))e2^«^ 



(38) 



Its formal solution is 



C{k,t) 



dt 



, UJk{t') 



2ujk 



Jy(l + 2AA+(£0) X 



(39) 



Substituting it into Eq. (p7| ), we obtain a closed equation 
for M+{k, t) similar to ||]: 



U+{k,t) = 



dt''^^(l + 2N+{k,t')] X 



2t^/c J-oo ~ ^kit') 

cos[2ef (t) - 2ef{t')] 



(40) 



Eq. ( [iof ) is a quantum kinetic equation which determines 
the time evolution of the number of particles of a fixed 
momentum k^ > k^. Note that the background field does 
not contribute to the kinetic equation directly, but only 
via the frequency of the quantum fluctuations. Therefore, 
the change of M+{k,t) in time in Eq. (|o|) is due to 
particle production during the fluctuations. 

In the regime of the negative frequency squared, when 

k^ < k1 and ujk = ±«t^fe = ±i\Jk1 — A:^, one of the phase 
factors in the ansatz (|^), T^{t) or Vt{t), grows exponen- 
tially in time. Instead of oscillations we have an expo- 
nential growth of long wavelength quantum fluctuations 
with momenta J? < fc^. This is the so-called tachyonic 
instabihty |l|-||. 

Such a tachyonic regime is realized for potential 
paramters a/fi^ < 1. Whether the system evolves in the 
tachyonic or non-tachyonic regime is dynamically fixed 
by the time dependent critical momentum: 



1 ^ I 



cos(0//) + a < 
otherwise 



(41) 



plotted in Fig. |i| for different parameters a//Lt^. For 
a/iJ? > 1 the critical momentum is zero since the fre- 
quency is always positive and no tachyonic modes can 
establish. The critical momentum for a//i^ < 1 oscil- 
lates in tune with the time dependence of the vacuum 
mean field (j). The time evolution shows that the same 
momentum state can change its nature during the evo- 
lution. In that case a different kinetic equation must be 
derived and solved which evolves all tachyonic modes in 
time. Therefore the analytical and numerical treatment 
is a complicated challenge and a quantitative analysis of 
a//i^ < 1 states will be provided elsewhere. 
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FIG. 1. The dependence of the critical momentum 
Pc = kc/fi, Eq. (W^), on time t — t /i. A non-vanishing 
value indicates the appearance of tachyonic modes. The time 
dependence of pc is due to the alternating (ji-field and depends 
strongly on the choice of the initial values. 

The total number of particles of all modes at any time 
t is given by 



(42) 



Simple power counting yields that this expression is fi- 
nite. 



V. DECAY OF THE CP ODD PHASE 

In the vicinity of the phase transition the QCD vacuum 
rearranges: chiral symmetry breaking and confinement 
drive quark matter into hadronic states. In this region 
topological phenomena such as the appearance of a non- 
vanishing QCD 6 angle may occur. As a consequence the 
CP symmetry is dynamically broken and in these regions 
a CP-odd phase can occur (paramater Set IV in Table 
I). In the rapid cooling of the hot QCD matter down 
to the critical temperature Td, the Uyi(l) breaking gets 
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restored either completely (Set I) or partially (Sets II, 
III). After such a quench of the effective potential (43), 
the system is in the false vacuum state. The decay of 
this CP odd phase is a time dependent process and can 
be study within the kinetic approach introduced in the 
previous section. 

The potential of the effective Lagrangian density (0) 



V{v/f) 



(43) 



is plotted in Fig. || for different values of the potential 
parameter, a//i^, accor ding to Table |. 




FIG. 2. The shape of the potential V{ri/f), is plotted 
for different values of a/n^. Local minima characteristic for 
a/fi^ < 1 assumed at large temperatures disappear due to an 
applied fast quench. 



Set 


I 


II 


III 


IV 


aj yi^ 1 (aj [P')^ac. 


1 


1/2 


1/4 






TABLE I. Different values oia/ as used in the numerical 
calculation, {a/ ^^)vac ~ 4.24 is the vacuum value for which 
all mesons have their vacuum masses. Set I assumes a fast 
quench after which the vacuum value is immediately reached, 
i.e. the r)' mass assumes its vacuum value in the vicinity of Td- 
This scenario is compared with slow quenches corresponding 
to parameters given in Set II and Set III. Set IV leads to 
the appearance of tachyonic modes; a value only possible for 
T > Td. 



Starting from the quark gluon plasma phase in which 
a/y? is suppressed. Set IV in Table |, the potential 
changes from the cosine shape to a parabolic shape due 
a sudden quench at the deconfinement phase transition. 
The metastable states located in the local minima of the 
potential roll smoothly back into the trivial minimum 
and oscillate around it. This situation is formalized in 
assumption ([ll|): The r/-field can be decomposed into its 
vacuum mean value 4>(t) and its quantum fluctuations x- 
During the decay, energy is transferred from (j) to x- 
a result, (j) is damped, while the number of particles in 
quantum fluctuations increases. It is assumed that dur- 
ing this process the temperature does not change essen- 
tially, and particle production proceeds in a fixed poten- 
tial characterized by a//i^. This process takes place on a 
time scale typical for the hadronisation process: 1 — 10 
fm/c. 




FIG. 3. The solution of the vacuum mean field equation as 
function of time, Eq. (Bl), is shown for different values of 
a/ {c.f. Table M) for the initial conditions 0(0)// = 27r and 
4>{0)/f — 0. Note that for a/fji^ = 0, <^(t) would be constant. 

The potential parameter a//Lt^ depends on temperature 
due to medium dependent meson masses. However its ex- 
act behaviour near the critical temperature is unknown. 
Lattice calculations as well as QCD models suggest that 
TT, K and rj meson properties have only a weak depen- 
dence on T. About the response of 77' to increasing T is 
much less known and therefore we explore different sce- 
narios summarised in Table |[ Set I assumes that the 
medium dependence is negligible. Set II (III, IV) corre- 
sponds to an in-medium reduction of the rj' mass of about 
20 (40, 60)% applying the simple equations given in con- 
nection with Eq. (^, 1^^. It is important to note that 
/ ~ /it can be considered as an order parameter for the 
chiral phase transition and hence is strongly suppressed 
at Td, i.e /(T - Ta) = 0.1 /(T = 0). In the herein 
applied scenario, the in-medium dependent parameters 
a/fj? and / change only during the fast quench. Their 
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time-dependence can therefore be neglected. 

The solution of the non- linear Eq. (Bl) for 0(r) with 
different a//i^ is plotted in Fig. ||, employing the initial 
conditions (/)(0)// = 2tt and {d(p{T) /dT)T-=o — through- 
out the numerical calculations. We see that the period 
and the amplitude of the oscillations of (/)(t) vary with 
the change of a/fi^. The oscillations are not damped 
since feedback of the fluctuations on the mean field is 
neglected. The field cj) provides the background field for 
the solution of the quantum kinetic equation given in Eq. 

(PD. 



frequencies, Eq. (B4), corresponding to a//i^ > 1 given 
in Table | 



VI. NUMERICAL RESULTS 

The solution of the quantum kinetic equation de- 
scribes the production of rj' particles: the momentum 
dependence and it's time evolution. The strong back 
ground field, (fi, leading to a sizeable particle production 
rate is given by the solution of the non- linear Eq. (15), 
see Fig. |[ 




Time 5^^~T>^ 2 

4 Momentum 

FIG. 4. The time evolution of the momentum distribu- 
tion function for parameter Set III. Most of the mesons 
are produced with small momenta but additional resonance 
bands appear for larger momenta; their maximal amplitude 
is smaller. The time evolution is characterized by an increase 
of the particle number and a repeated spike structure. 

We perform the numerical calculation using dimension- 
less variables and solve the kinetic equat ion as a system 
of coupled differential equations, Eqs. (B8-B1C), intro- 
duced in Appendix B. The decay starts at r = 0, for 
which Af+{p, 0) = 0; 0(0) = 2nf and 0(0) = define the 
initial conditions. 

As result we obtain the number of particles produced 
during the decay of the CP odd phase in the false vacuum. 
Herein we restricted ourselves to the study of the non- 
tachyonic regime, i.e. we explore the solution for positive 
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FIG. 5. The time evolution of the particle number for two 
different aju^ > 1 when the system is in the non-tachyonic 
regime, ( [B3|) . The double spike structure on top of the rapid 
growth repeats periodically in tune with the mean field's fre- 
quency (c./. Fig. 1^. An estimate in low density approxima- 
tion shows that inclusion of the Bose quantum statistics leads 
to a pronounced enhancement. 

In Fig. I we show the complete numerical solution. 
Two features are apparent: (i) the fast increase is char- 
acterised by a repeated structure on top of the curve, (ii) 
additional to the occupation of low momentum states we 
observe the appearance of resonance bands at larger mo- 
menta. 

In Fig. H we plot the time evolution of the particle 
number for zero momentum and compare the solution for 
Set I and Set III. We observe a very fast increase of the 
number of produced particles. A maximum occupation 
of a given momentum state at a given time is reached 
for small values of the potential parameter, i.e. using Set 
III. For larger values of the potential parameter, e.g. Set 
I, less particles are produced in a given time since the 
source term is suppressed by a larger mass term a//i^ in 
Lu{p), see (|B4|). 



The most striking feature in this plot is the periodically 
repeated spike structure on top of the overall growth, {c.f. 
p^). This pattern appears with the same frequency as 
the background field oscillates, see Fig. ^. When back 
reactions are included this would possibly not be the case. 
The spike structure is smoother for Set I compared to Set 
III but still characteristic for the evolution. 

Herein we also compare the full solution with the low 
density approximation (l.d.). The low density approxi- 
mation assumes that J\f+ {p, t) <C 1 and hence suggests 
that the solution of the kinetic equation does not depend 
on the pre-history of the systems evolution. Any calcu- 
lation which does not retain quantum statistical effects 
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necessarily employs this ansatz. From Fig. ^ it is plain 
that the inclusion of the Bose statistical factor into the 
kinetic equation leads to Bose enhancement as soon as 
~ 1, appearing at t 1. This effect becomes more 
pronounced with increasing time. 
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FIG. 6. The particle number as function of momentum, 
p = \p\, for two different a/n^ > 1. Bose enhancement of 
mainly the low momentum states is apparent. A characteris- 
tic second resonance band appears for large momenta. 

The momentum dependence at a given time, Fig. 
shows that most of the particles are produced with small 
momenta. Additional resonance bands appear. The 
smaller the value of the potential parameter is reached in 
the quench the closer the second maxima appears to the 
first one. The reason for this resonance effect is typical 
for the Mathieu type equation: the two intrinsic frequen- 
cies of the background field and of the production process 
are of the same order of magnitude and resonances are 
likely to appear. 

In Fig. ^ we also compare the momentum dependence 
of the full non-Markovian solution with a calculation in 
low density limit. It is apparent that the Bose enhance- 
ment acts naturally on the lower momenta. For the con- 
sidered case the occupation number is enhanced by a 
factor of 4. For large momenta details of the quantum 
statistics are suppressed. Therefore the higher resonance 
bands are much less affected by quantum corrections. It 
is plain from this study that quantum statistical effects 
cannot be neglected: the low density approximation is 
invalid if the produced number density exceeds a critical 
value at very early times of the evolution. 



VII. SUMMARY 

Starting from the singlet Witten-DiVecchia-Veneziano 
effective Lagrangian we have derived a quantum kinetic 
equation describing the production of yy'-mesons from a 



CP-odd metastable vacuum state. We have employed a 
general method based on the evolution operator holding 
also for other model langrangians. The vacuum mean 
field provides a classical, self-interacting strong back- 
ground field. Quantum fiuctuations around the dynam- 
ical mean field value are considered but their feedback 
to the background field is neglected. Due to these quan- 
tum fluctuations particles are produced and the time evo- 
lution of this process is described by a non-Markovian 
equation for the distribution function of the produced rj' 
mesons. 

We find that the details of the decay process depend 
strongly on the applied quench. The number of pro- 
duced particles is much larger when the rf mass is sup- 
pressed in the vicinity of the phase boundary. Most of 
the particles are produced with low momenta, for large 
momenta additional resonances appear. Furthermore, 
we have demonstrated that quantum statistical effects 
are important and lead to a pronounced enhancement of 
the particle occupation number for low momenta. In the 
case a//i^ < 1, tachyonic instabilities occur for momenta 
smaller than a critical value. This regime has not been 
considered herein. 

The numerical investigation of the tachyonic modes 
and the inclusion of back reactions promise further in- 
sight into the decay of CP odd metastable states and its 
realization will be reported elsewhere. 
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APPENDIX A: IN-FIELD QUANTIZATION 

The in-field is a solution of the equation 

(□ -f- ml)r],n = 0. (Al) 

The in-field operators fulfill periodic boundary conditions 
and are expanded in Fourier modes 

k 

7r,„(f,t) = -^^e-''=^7r„(fc,t), (A3) 

k 

where the summation is over discrete momenta k — '^n, 
(ni,n2,n3) and 
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r?„(fc,t) = rO-(t)a„(fc) + r°'*(t)aL(-fc). 



(A4) 
(A5) 



r'l{t)aU-k)-T':'*{t)al{k) 



The time-independent creation and annihilation opera- 
tors obey the commutation relations 



[a„(fc),aL(fc')]-='5,,t„ 



(A6) 



all other commutators vanish. The function ^^{t) is given 
by 



xp{-«t^fci} 



(A7) 



with uj^ = \J + rriQ. Since the field r]{x,t) is real, we 

have ?7,L(fc,i) = T]in{-k,t), and 7r|„(fc, i) = 7rj„(-fc,i). 
The vacuum state |0; in) = |0) is defined as vanishing 
under the action of the annihilation operators ai„(fc)|0) = 
0. 



Eq. ( B3 ) is an integro-differential equation. It can be 
re-expressed by introducing 

u{p, r) EE r dr'^ir') (l + 2U+{p, r')) x 
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sin[2ef (r) - 2Qf{T') 



(B5) 

v{p, r)^ i dr'^ir') (l + 2Af+{p, r')) x (B6) 
Jo ^ ' 

coA2%\t)-2%\t'% (B7) 



with the initial conditions 0) — 0) = 0, in which 
case we have 



(B8) 



d_ 
dr 
d_ 
d7 



v{p, t) = ^(i + 2M+{p. t)) - 2(ip«(p, r), (B9) 



u{p,t) = 2ujpv{p,T). 



(BIO) 



APPENDIX B: NUMERICAL REALIZATION 



The evolution of in the decay is governed by equa- 
tion (p^). Introducing the dimensionless vacuum mean 
field (j)/ f and the dimensionless time variable r = /it, we 
rewrite (Esl) as 



= 0, 



(Bl) 



with the one parameter a/ p^ characterizing the s olut ion. 

Note that for small a/n^ « one can replace (Bl) by 
the Sine- Gordon equation 



dr2 V 



Mr) 
f 



Mr) 
f 



0, 



(B2) 



the subscript (0) in (/)(t) indicates the zero value of a/fi^. 
The solution of Eq. (B2) is a Jacobian elliptic function. 
Herein we do not make this approximation and solve (Bl) 
numerically for nonzero values of a/ fi^. 

For the numerical study we introduce dimensionless 
variables for the kinetic equations and obtain 



-^AA+(p, T)^^ I ■ dr'^^ir') ( 1 + 2AA+(p, r'] 



dr 



20j. 



P Jo 



cos[2e;\r) - 2e;'{r')] 
where the dimensionless frequency is 
1 



P 



f 



(- 



(B3) 



(B4) 



with p^ = (P//i2). 
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